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Abstract 

In this paper we describe an algorithm that takes as input a description of a semi- 
algebraic set S C R'^, defined by a Boolean formula with atoms of the form P > 
0,P < 0,P = for P G P C R[Xi, . . . and outputs the first f + 1 Betti 

numbers of S, bQ{S), . . . , be{S). The complexity of the algorithm is (sd)'^°**\ where 
where s = #(7-") and d = maxpgp deg(P), which is singly exponential in k for 
i any fixed constant. Previously, singly exponential time algorithms were known 
only for computing the Euler-Poincare characteristic, the zero-th and the first Betti 
numbers. 
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1 Introduction 



Let R be a real closed field and S" C R'^ a semi-algebraic set defined by a 
Boolean formula with atoms of the form P > 0,P < 0,P = for P G 
V C R[Xi, . . . , Xk] (we call such a set a P-semi- algebraic set and the corre- 
sponding forrnula a P -formula). It is well known ()0]einik5lt l0P49l : lMi1nor64l : 



sponamg lorrnula a /- ^-lormula) . it is weii Jinown i|l JiemiKo it IUr4Ht IMnnorr)4l ; 
Thorn 6,4 iRoi lGVn5h that the topological complexity of S (measured by the 



various Betti numbers of 5*) is bounded by {sd)'^^''\ where s = if{V) and 
d = maxpgp deg(P). Note that these bounds are singly exponenti al in k . More 
precise bounds on the individual Betti numbers of S appear in (|B03I ). Even 
though the Betti numbers of S are bounded singly exponentially in k, there 
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is no known algorithm for producing a singly exponential sized triangulation 
of S (which would immediately imply a singly exponential algorithm for com- 
puting the Betti numbers of S). In fact, designing a singly exponential time 
algorithm for computing the Betti numbers of semi-algebraic sets is one of 
the outstanding open problems in algorithmic semi-algebraic geometry. More 
recently, determining the exact complexity of computing the Betti numbers 
of semi-algeb raic set s has attracted the attention of computational complex- 
ity theorists (jBC04^ . who are interested in developing a theory of counting 
complexity classes for the Blum-Shub-Smale model of real Turing machines. 



Doubly exponential algorithms (with complexity (s(i)^°**') for computing all 
the Betti numbers are known, since it is possible to obtain a triangulation 
of S in doubly ex ponential time using cylindrical algebraic decomposition 
f)Collins7,4lBPR,n3t ). In the absence of singly exponential time algorithms for 
computing triangulations of semi-algebraic sets, algorithms with single ex- 
ponential complexity are known only for the problems of testing emptiness 
f)Renegar92l : lBPR9.5l ). computing the zero-th Betti numbe r (i.e. the number 
of semi- algebraically connected components of S) (IGV92I: Cannv93l : GR92I : 
BPR93), as well as the Euler-Poincare characteristic of S (|B99|). Very recently 



a singly exponential time algorithm has been developed for t he probl em of 
computing the first Betti number of a given semi-algebraic set ()BPR04| ) . 



In this paper we describe an algorithm, which given a family V C R[Xi, . . . , X^], 
a P- formula describing a P-semi-algebraic set S* C R^, and a number i,0 < 
i < k as input, outputs the first i Betti numbers of S. For constant £, the 
complexity of the algorithm is singly exponential in k. We remark that using 
Alexander duality, we immediately get a singly exponential algorithm for com- 
puting the top i Betti numbers too. However, the complexity of our algorithm 
becomes doubly exponential if we want to compute the middle Betti numbers 
of a semi- algebraic set using it. 



There are two main ingredients in our algorithm for computing the first i 
Betti numbers of a given closed semi- algebraic set. The first ingredient is a 
result proved in (jBPRO^, which enables us to compute a singly exponential 
sized cover of the given semi-algebraic set consisting of closed, contractible 
semi-algebraic sets, in single exponential time. The number and the degrees of 
the polynomials used to define the sets in this cover are also bounded singly 
exponentially. 



The second ingredient, which is the main contribution of this paper, is an 
algorithm which uses the covering algorithm recursively and computes in singly 
exponential time a complex whose cohomology groups are isomorphic to the 
first i cohomology groups of the input set. This complex is of singly exponential 
size. 
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The main result of the paper is the following. 

Main Result: For any given £, there is an algorithm that takes as input a V- 
formula describing a semi-algebraic set S C R'^, and outputs hQ^S), . . . , bi{S). 
The complexity of the algorithm is {sdY°^'\ where s = and d = 

maxpg-p deg(P). Note that the complexity is singly exponential in k for every 
fixed L 

The paper is organized as follows. In Section 2, we recall some basic defi- 
nitions from algebraic topology and fix notations. In Section 4 we describe 
the construction of the complexes which allow us to compute the the first ^ 
Betti numbers of a given semi-algebraic set. In Section 5 we recall t he inputs , 



outputs and complexities of a few algorithms described in detail in (jBPROJ), 
which we use in our algorithm. In Section 6 we describe our algorithm for 
computing the first £ Betti numbers, prove its correctness as well as the com- 
plexity bounds. Finally in Section 7 we comment on issues related to practical 
implementation. 



2 Mathematical Preliminaries 



In this section, we recall some basic facts about semi-algebraic sets as well as 
the definitions of complexes and double complexes of vector spaces, and fix 
some notations. 



2.1 Semi- algebraic sets and their cohomology groups 



Let R be a real closed field. If P is a finite subset of R[Xi, . . . , X^], we write 
the set of zeros of V in R^ as 

Z(P, R^) = {xeR^\ f\ P{x) = 0}. 

Per 

We denote by B{0,r) the open ball with center and radius r. 

Let Q and V be finite subsets of R[Xi, . . . , X^], Z = Z(Q, R^), and Zr = 
Zr\B{0, r). A sign condition on V is an element of {0, 1, — l}^- The realization 
of the sign condition a over Z, TZ{(t, Z), is the basic semi-algebraic set 

{x G R'^ I A Qi^) = /\ A sign(P(x)) = a{P)}. 

Q£Q PeV 

The realization of the sign condition a over Zr., TZ{a,Zr), is the basic semi- 
algebraic set 7?.((T, Z) n B{0, r). For the rest of the paper, we fix an open ball 
5(0, r) with center and radius r big enough so that, for every sign condition 
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cr, 7^(cr, Z) and 7l{a, Zr) are homeomorphic. This is al ways p ossible by the 
local conical structure at infinity of semi-algebraic sets f ()BCR[ l. page 225). 



A closed a nd bound ed semi-algebraic set S" C R'^ is semi-algebraically triangu 



lable (see flBPR03[ )). and we denote by H*(5') the i-th simplicial cohomology 



group of S with rational coefficients. The groups H*(S') are invariant under 
semi-algebraic homeomorphisms and coincide with the corresponding singular 
cohomology groups when R = M. We denote by bi{S) the i-th Betti number 
of 5* (that is, the dimension of H*(5') as a vector space), and b{S) the sum 
J2i bi{S). For a closed but not necessarily bounded semi-algebraic set S C R'^, 
we will denote by H*(5') the i-th simplicial cohomology group of S" fl B{0, r), 
where r is sufficiently large. The sets Sr\B{0, r) are semi-algebraically homeo- 
morphic for all sufficiently large r > 0, by the local conical structure at infinity 
of semi- algebraic sets, and hence this definition makes sense. 

The definition of cohomology groups of arbitrary semi-algebraic sets in R'^ 
requires some care and several possibilities exist. In this paper, we follow 
()BPR,03[ ) and define the cohomology groups of realizations of sign conditions 
as follows. 

Let R denote a real closed field and R' a real closed field containing R. Given 
a semi-algebraic set S in R'^, the extension of S to R', denoted Ext(S', R'), is 
the semi-algebraic subset of R''^ defined by the same quantifier free formula 
that defines S. The set Ext(5', R') is well defined (i.e. it only depends on the 
set S and not on the quantifier free formul a chosen to describe it). This is an 
easy consequence of the transfer principle 



Now, let 5 C R^ be a P-semialgebraic set, where V = {Pi,...,Ps} is a 
finite subset of R[Xi, . . . ,Xk]. Let (f){X) be a quantifier-free formula defining 
S. Let Pi = So ai,Q,X° where the ai^a € R- Let A = {. . . , Ai^a, ■ ■ ■) denote 
the vector of variables corresponding to the coefficients of the polynomials 
in the family V, and let a = (. . . , ai^a, • • •) ^ R^ denote the vector of the 
actual coefficients of the polynomials in V. Let ip{A,X) denote the formula 
obtained from 4>{X) by replacing each coefficient of each polynomial in V by 
the corresponding variable, so that 4>{X) = ijjia ^ X). It fo llows from Hardt's 
triviality theorem for semi-algebraic mappings ( HardtSol ). that there exists. 



a' G M^g such that denoting by S' C M^ig the semi- algebraic set defined by 
'?/'(a', X), the semi-algebraic set Ext(S", R) has the same homeomorphism type 
as S . Here, M^ig is the field of real algebraic numbers. We define the cohomology 
groups of S to be the singular cohomology groups of Ext (5", R). It follows from 
the Tarski-Seidenberg transfer principle, and the corresponding property of 
singular cohomology groups, that the cohomology groups defined this way are 
invariant under semi-algebraic homotopies. It is also clear that this definition is 
compatible with the simplicial cohomology for closed, bounded semi-algebraic 
sets, and the singular cohomology groups when the ground field is R. Finally 
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it is also clear that, the Betti numbers are not changed after extension: 

k{S)=k{Ext{S,R')). 



Note that we define the co-homology groups of arbitrary semi-algebraic sets 
as above in order to treat semi-algebraic sets over arbitrary (possibly non- 
archimedean) real closed fields R, for which the standard proofs of the ho- 
mology axioms (in pa rticular the e xcision axiom) break down for singular 
homology groups (see (|Knebusch89l ). page XIII). If one is only interested in 
the case, R = M, then singular co-homology groups suffice. 



2.2 Complex of Vector Spaces 



A sequence {C^}, p E Z, of Q-vector spaces together with a sequence {S^} of 
homomorphisms 6^ : —>■ C^"*"^ (called differentials) for which 5^"^ 5^ = 
for all p is called a complex. When it is clear from context, we will drop the 
supercripts from the differentials for the sake of readability. 

The cohomology groups, W{C') are defined by, 

W{C') = ZP{C*)/BP{C*), 

where 5^(0*) = lm{6P-^), and Zp{C') = Ker(5P) and we will denote by 
H*(C') the graded vector space 0pHP(C'). 

The cohomology groups, W{C*), are all Q-vector spaces (finite dimensional if 
the vector spaces C^'s are themselves finite dimensional). We will henceforth 
omit reference to the field of coefficients Q which is fixed throughout the rest 
of the paper. 



2.3 Homomorphisms of Complexes 



Given two complexes, C* = (C^,5^) and D* = (D^,^^), a homomorphism of 
complexes, 0* : C* ^ D*, is a sequence of homomorphisms (p^ : C'p —>■ for 
which 6P (jf = (ff+^ SP for all p. 



In other words, the following diagram is commutative. 



CP 
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A homomorphism of complexes, 0* : C* — > D*, induces homorphisms, 0* : 
H*(C*) H'(D*) and wc will denote the corresponding homomorphism be- 
tween the graded vector spaces H*(C*), H*(D*) by (p*. The homomorphism 0* 
is called a quasi-isomorphism if the homomorphism 0* is an isomorphism. 

Given two complexes C* and D*, their direct sum denoted by C*©D*, is again a 
complex with itsp-th term being C^©D^. Moreover, given two homomorphisms 
of complexes, 

0* : C ^ C*, 
V'' : D* ^ D', 

their direct sum 

(/)• © V : C* e ^ C © D", 

is again a homomorphism of complexes defined componentwise. Note that if 
we specify a basis for the different terms of the complexes C*,C*,D*,D*, as 
well as the matrices for the homomorphisms (p*, ip* then we can write down the 
matrix for the direct sum homomorphism 0* © -0' as a sum of block-matrices 
using elementary linear algebra. 

2.4 The Nerve Lemma and Generalizations 

We first define formally the notion of a cover of a closed, bounded semi- 
algebraic set. 

Definition 2.1 Let S G K'' be a closed and bounded semi-algebraic set. A 
cover, C{S), of S consists of an ordered index set, which by a slight abuse of 
language we also denote by C{S), and a map that associates to each a G C{S), 
a closed and bounded semi-algebraic subset Sa C S, such that S — Lia£C(s)Sa- 

For do, . . . , Op, G C{S), we associate to the formal product, ckq • • • CKp, the closed 
and bounded semi-algebraic set Sao-ap = Sao H ■ • • fl Sa^- 

Recall that the 0-th simplicial cohomology group of a closed and bounded 
semi-algebraic set X, H''(X), can be identified with the Q- vector space of Q- 
valued locally constant functions on X. Clearly, the dimension of H'^(X) is 
equal to the number of connected components of X. 

For cco, CKi, . . . , CKp, (3 G C{S), and /? {ckq, . . . , ctp}, let 

be the homomorphism defined as follows. Given a locally constant function, 
(p G H°(5'ao-..aj,), raa-ap;0{,(p) IS the locally constant function on Sao-ap-p ob- 
tained by restricting (p to Sao-ap-p- 
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We define the generalized restriction homomorphisms, 

ao<—<ap,aieC{S) ao<—<ap+i,aieC{S) 

by 

<^^(0)q;o---q;p+i 5Z ("■'■) '"Q;o---cEr--o;p+i;ai(0Q;o---cfi---ap+i)' (■'■) 

0<i<p+l 

where e 0ao<-<apGC(S) H°(5'„o...aJ and 

''"ao---cfi---ap+i;ai -^^ ^■'^^ restriction ho- 
momorphism defined previously. The sequence of homomorphisms 5^ gives rise 
to a complex, L*(C(5')), defined by, 

L"(C(^)) = H°(5„„...„J, 

ao<---<ap,ai&C{S) 

with the differentials 5^ : Lp{C{S)) U+\C{S)) defined in (1). The complex 
L'{C{S)) is often referred to as the nerve complex oi the cover C{S). 

For any £ > 0, we will denote by L*{C{S)) the truncated complex, defined by. 



mcis))^LP{c{s)),o<p<e, 

0, p>t 

Notice that once we have a cover of S, and we identify the connected compo- 
nents of the various intersections, Sao-ap, we have natural bases for the vector 
spaces 

ao<---<ap,aieC{S) 

appearing as terms of the nerve complex. Moreover, the matrices correspond- 
ing to the homomorphisms in this basis, depend only on the inclusion rela- 
tionships between the connected components of Sao-ap+i and those of Sao-ap- 

We say that the cover C{S) satisfies the Leray property if each non-empty 
intersection Sao-ap is contractible. Clearly, in this case 

H°(-S„) =Q, if5ao-ap7^0 
^ 0, if Sao-ap = 0. 

It is a classical fact (usually referred to as the nerve lemma) that. 

Theorem 2.2 (Nerve Lemma) Suppose that the cover C{S) satisfies the Leray 
property. Then for each i > 0, 

ff (L'(C(>5))) ^ R'iS). 
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Proof: See (|Rotmanl ) . 



□ 



Thus, Theorem 2.2 gives a method for computing the Betti numbers of 5* using 
hnear algebra, from a cover of S by contractible sets for which all non-empty 
intersections are also contractible, once we are able to test emptiness of the 



various intersections Sao-.-ap- 



Now suppose that each individual member, Sa^ of the cover is contractible, 
but the various intersections Sao-ap are not necessarily contractible for p > 1. 
Theorem 2.2 does not hold in this case. However, the following is proved in 



(jBPROJ) 



Theorem 2.3 Suppose that each individual member, Sag of the cover C{S) is 
contractible. Then, 

ff (L;(C(5))) - R^S), 

for i = 0,1. 



Proof: See (|BPROJ). □ 



Notice that from a cover by contractible sets. Theorem 2.3 allows us to com- 
pute using linear algebra, bo^S) and bi{S), once we have identified the non- 
empty connected components of the pair-wise and triple-wise intersections of 
the sets in the cover, and their inclusion relationships. It is quite easy to see 
that if we extend the complex in Theorem 2.3 by one more term, that is con- 
sider the complex, L3(C(5')), then the cohomology of the complex does not 
yield information about H^(S'). Just consider the cover of the standard sphere 
52 c R^ and the cover {Hi,H2} of S"^ where Hi,H2 are closed hemispheres 
meeting at the equator. The corresponding complex, L'{C), is as follows. 



^ H°(ifi) H°(i72) ^ H°(/fi n iJa) ^ ^ 



Clearly, H^(L3(C(S'))) 9^ H^(S'^), and indeed it is impossible to compute bi{S) 
just from the information on the number of connected components of inter- 
sections of the sets of a cover by contractible sets for, i > 2. For example, the 
nerve complex coresponding to the cover of the sphere by two hemispheres is 
ismorphic to the nerve complex of a cover of the unit segment [0, 1] by the 
subsets [0, 1/2] and [1/2, 1], but clearly li^iS^) = Q, while B.^{[0, 1]) = 0. 

In order to deal with covers not satisfying the Leray property, it is necessary 
to consider a generalization of the nerve complex, namely a double complex 
arising from the generalized Mayer- Vietoris exact sequence. The construction 
of this double complex (which is quite classical) in fact motivates the design 
of our algorithm, which we describe in detail in Section 6. 
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3 Mayer- Vietoris 



3.1 Double Complexes 

In this section, we recall the basic notions of a double complex of vector 
spaces and associated spectral sequences. A double complex is a bi-graded 
vector space, 

with co-boundary operators d : C^'^ and 5 : C'*'''' Q,p+^^i and such 

that d5 -\- 5d — We say that C*'* is a first quadrant double complex, if it 
satisfies the condition that C^'^ = if either p < or ^ < 0. Double complexes 
lying in other quadrants are defined in an analogous manner. 

The complex defined by 

Tot" (€•'•) = C^''«, 

p+q=n 

with differential 

D" = d± 5 : Tot"(C''*) — > Tot"+^(C*''), 
is denoted by Tot'(C*'*) and called the associated total complex ofC*'*. 

3.2 Spectral Sequences 

A spectral sequence is a sequence of bigraded complexes {Ej.,dr : E'^'"^ 
£jp+r,q-r+r'j ^^^^ Figurc 1) such that the complex Er+i is obtained from E,. by 
taking its cohomology with respect to dr (that is Er+i = B.d^{Er)). 

There are two spectral sequences, "E'^''', (corresponding to taking row- 
wise or column-wise filtrations respectively) associated with a first quadrant 
double complex C*'', which will be important for us. Both of these converge to 
H*(Tot'(C''*)). This means that the homomorphisms, dr are eventually zero, 
and hence the spectral sequences stabilize, and 

'E'J - "E'^' = ff (Tof(C-'-)), 

p+q=i p+q=i 

for each i > 0. 

The first terms of these are 

'£;i = H,(C-'-),% = H,H5(C-'-), 
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and 



p+q=l p+q=e+l 

Fig. l.dr-. ^ 



Given two (first quadrant) double complexes, C*'* and C*'*, a homomorphism 
of double complexes, 



is a collection of homomorphisms, cff''^ : C^'^ — C^'^, such that the following 
diagrams commute. 



QP,<1 1 QP-^+l 



(7P,? ° ) (jp,g+i 
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A homomorphism of double complexes, 

: C'' > C*'*, 

induces an homomorphism of the corresponding total complexes which we will 
denote by, 

Tof(0*'') : Tof(C''') — > Tof(C*''). 

It also induces homomorphisms, '0s • 'Es — > 'Eg (respectively, "(pg '■ "Eg — > 
"Eg) between the associated spectral sequences (corresponding either to the 
row-wise or column-wise filtrations) . For the p recise definition of homomor- 
phisms of spect ral sequenc es, see ()MclearvOll ). We will need the following 
useful fact (see flMclearvOlh . page 66, Theorem 3.4 for a proof). 

Proposition 3.1 //Vs (respectively, "(pg) is an isomorphism for some s > 1, 
then 'E^j.''^ and'E^''^ (repectively, "E^''^ and"El''' ) are isomorphic for allr > s. 
In particular, the induced homomorphism, 

Tof(0''*) : Tof(C*'*) — ^ Tof(C"'*) 

is a quasi-isomorphism. 



3.3 The Mayer-Vietoris Double Complex 



Let Ai,. . . ,An be sub-complexes of a finite simplicial complex A such that 
A = Ai U ■ ■ ■ U An- Note that the intersections of any number of the sub- 
complexes, Ai, is again a sub-complex of A. We will denote by Aa^^.-ap the 
sub-complex A^g fl ■ ■ ■ fl Aa^- 

Let C*(A) denote the Q- vector space of i co-chains of A, and C*{A), the 
complex 

> C'-'iA) ^ ^(A) ^ C^+\A) 

where d : C(yl) —>■ C'^^^(A) are the usual co-boundary homomorphisms. More 
precisely, given u G C{A), and a g + 1 simplex [ao, . . . , flg+i] G A, 

duj{[ao,. . .,ag+i]) = J2 (-l)'w([ao, • • • , di, . . . , Og+i]) (2) 

0<i<q+l 

(here and everywhere else in the paper " denotes omission). Now extend duj to 
a linear form on all of Cq+i{A) by linearity, to obtain an element of C'^'^^{A). 

The connecting homomorphisms are "generalized" restrictions and are defined 
below. 

The generalized Mayer-Vietoris sequence is the following exact sequence of 
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vector spaces. 



l<ao<n l<ao<a\<n 
l<ao<---<ap<n l<ao<---<ap+i<n 

where r' is induced by restriction and the connecting homomorphisms 5^'' are 
as follows. 

Given an u; G ®ao<---<ap C'^{^ao---ap) we define 6^''^{uj) as follows: 

First note that 6P''^uj G ®ao< --<ap+i i^ao -ap+i) , and it suffices to define 

for each (p + 2)-tuple 1 < < ■ ■ ■ < ctp+i < Note that, (5^'^ci;)ao,...,ap+i 
is a linear form on the vector space, Cg{Aaf,...a ^J, and hence is determined 
by its values on the g-simplices in the complex Aag...ap^-^- Furthermore, each 
g-simplex, s G Acn^...a^^^ is automatically a simplex of the complexes 



Aao-a,-ap+i, < 2 < p + 1. 



We define, 



ao,...,aj,...,ap+i 



o<i<p+i 



The fact t hat the g eneralized Mayer- Vietoris sequence is exact is classical (see 



( RotmarJ l or (|B03l ) for example). 



We now define the Mayer- Vietoris double complex of the complex A with 
respect to the subcomplexes Aa^, ^ < ao < n, which we will denote by A/"*'* (A) 
(we suppress the dependence of the complex on sub-complexes Aa^ in the 
notation since this dependence will be clear from context). 

Definition 3.2 The Mayer- Vietoris double complex of a simplicial complex 
A with respect to the subcomplexes Aa^,! < ao < n, Af*'*{A), is the double 
complex defined by, 

= C^(A.o-.J- 

l<Q:o<--<ap<" 

The horizontal differentials are as defined above. The vertical differentials are 
those induced by the ones in the different complexes, C {A^Q.-ap) ■ 
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A^*'*(y4) is depicted in the following figure. 



I 



ao 



ao<ai 



ao<ai 



ao<ai 



For ani/ 1 > 0, tfe denote by A/"*'* (A) t/ie following truncated complex. 

Mr{A) = AfP^'^iA), 0<p + q<t, 
Xt%A) = 0, otherwise. 



The following proposition is classical (see (|Rotman[ ) or (|B03[ ) for a proof) and 
follows from the exactness of the generalized Mayer- Vietoris sequence. 

Proposition 3.3 The spectral sequences, ' E^," E^, associated to M*'* {A) con- 
verge to W{A) and thus, 

H*(Tof(Ar'''(A))) = B*{A). 

Moreover, the homomorphism 

^' : C'{A) Tof(Ar'''(A)) 

induced by the homomorphism r* (in the generalized Mayer- Vietoris sequence) 
is a quasi-isomorphism. 

We denote by C'_^]^(A) the truncation of the complex C*{A) after the (£+ l)-st 
term. As an immediate corollary we have that, 



Corollary 3.4 For any i > 0, the homomorphism 

^;+i:C,V(A)^Tof(Ar;;;(A)) 



(3) 



induced by the homomorphism r' (in the generalized Mayer- Vietoris sequence) 
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is a quasi-isomorphism. Hence, for < i < £, 

ff(Tof(Ar,-;;(yl))) = ff(A). 

Remark 3.5 Notice that in the truncated Mayer-Vietoris double complex, 
M''*{A), the 0-th column is a complex having at most t + 1 non-zero terms, the 
first column can have at most t non-zero terms, and in general the i-th column 
has at most t + 1 — i non-zero terms. This observation plays a crucial role in 
the inductive argument used later in the paper (in the proof of Proposition 
4.3). 



4 Double complexes eissociated to certain covers 

We begin with a definition. 

Definition 4.1 Let V he a finite subset o/R[Xi, . . . ,X]^. A V -closed formula 
is a formula constructed as follows: 

For each P &V, 

P = 0,P > 0,P < 0, 

are V-closed formulas. 
If $1 and $2 dfe V-closed formulas, $i A $2 and $1 V $2 are V-closed for- 
mulas. 

Clearly, TZ{^) — {x <Z \ ^{x)}, the realization of a V-closed formula 
is a closed semi- algebraic set and we call such a set a V-closed semi- algebraic 
set. 

In this section, we consider a fixed family of polynomials, V C Rf-'^i, . . . , ^fe], 
as well as a fixed P-closcd and bounded semi-algebraic set, ^S" C R*^. We also 
fix a number, i,0 < i < k. 

We define below (in Section 4.1) a finite set of indices, A5, which we call the set 
of admissible indices, and a map that associates to each a e A5 a closed and 

bounded semi-algebraic subset Xa C 5", which we call an admissible subset. 
To each a G A^, we associate its level, denoted level(Q;), which is an integer 
between and £. The set A5 will be partially ordered, and we denote by 
ancestors(Q;) C A5, the set of ancestors of a under this partial order. For 
a,P E A5, P e ancestors (a), implies that C Xp. 

For each admissible index a e Ag, we define a double complex, M.*'*{a), such 
that 

W{Tof{M'''{a))) ^ W{X^), 0<i<e- level(a). 
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The main idea behind the construction of the double complex Ai'''{a) is as 
follows. Associated to any cover of there exists a double complex (the 
Mayer- Vietoris dou ble co mplex) arising from the generalized Mayer- Vietoris 



exact sequence (see ( BOSh ). If the individual sets of the cover of X are all con- 



tractible, then the first column of the Mayer- Vietoris double complex is zero 
except at the first row. The cohomology groups of the associated total complex 
of the Mayer- Vietoris double complex are isomorphic to those of Xa and thus 
in order to compute bo^Xa), . . . ,bi_icvei{a){Xa), it suffices to compute a suit- 
able truncation of the Mayer- Vietoris double complex. However, computing 
(even the truncated) Mayer- Vietoris double complex directly within a singly 
exponential time complexity is not possible by any known method, since we 
are unable to compute triangulations of semi-algebraic sets in singly exponen- 
tial time. However, making use of the cover construction recursively, we are 
able to compute another double complex, Ai*'*{a), which has much smaller 
size but whose associated total complex is quasi-isomorphic to the truncated 
Mayer- Vietoris double complex and hence has isomorphic cohomology groups 
(see Proposition 4.6 below). The construction of A4*'*{a) is possible in singly 
exponential time since the covers can be computed in singly exponential time. 

Finally, given any closed and bounded semi- algebraic set X C R^, we wil l 



denote by C'{X), a fixed cover of X (we will use the construction in (jBPROJ) 
to compute such a cover). 



4-1 Admissible sets and Covers 



We now define As, and for each a G A5 a cover C{a) of X^ obtained by 
enlarging the cover C'{Xa)- 

Definition 4.2 (Admissible indices and covers) We define As by induction 
on level. 

(1) Firstly, G As, level(O) = 0, Xq = S, and C(0) = C'{S). The admissible 
indices at level 1 consists of all formal products, jd = ■ ai ■ ■ ■ aj-i ■ Oij, 
with ai G C(0) and < j < i + 1, and we define the associated semi- 
algebraic set by, 

Xf3 = Xao n ■ ■ ■ n Xaj ■ 

For each {ao, . . . , am} C {Po, . . . , /?„} C C{0), with n < £ + 1, 
ao ■ ■ ■ am £ ancestors(/5o ■ ■ ■ Pn), 

and G ancestors(/3o ■ ■ ■ Pn) ■ 

(2) We now inductively define the admissible indices at level i + 1, in terms 
of the admissible indices at level < i. For each a & As at level i, we 
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define C{a) as follows. Let ancestors(Q;) = {cci, . . . , a;^^}. Then, 

where [j denotes the disjoint union. All formal products, (3 — ao-ai - ■ -aj, 
with ai e C{a) and 0<j<£ — i + 1 are in As, and we define 

Xfj = n • • • n Xaj , 

and level(/3) — i + 

For each {ao, . . . , am} C {/3o, ■ ■ ■ , Pn} C C{a), with n < i — i + 1, 

tto ■ ■ ■ «m £ ancestors (/3o • • • Pn), 
and a G ancestors (/3o ■ ■ ■ Pn)- 

Moreover, for a' G C'{j3i /3jv ■ (^), each j3i is an ancestor of a' . We 

transitively close the ancestor relation, so that ancestor of an ancestor 
is also an ancestor. Moreover, if a^ - • • am, Po - ■ ■ Pn £ O'f^ such that 
for every j G {1, . . . , n} there exists i G {1, . . . , m} such that ai is an 
ancestor of j3j, then a^ - ■ ■ is an ancestor of Pq - ■ ■ (3^- 

Finally, the set of admissible indices at level i + 1 is 

U„eAs,level(a)=^^«0 | G C{a), < j < £ - i + 1} 

Observe that by the above definition, if a, /? G As and P G ancestors(a), then 
each a' G C{a) has a unique ancestor in each C{P), which we will denote by 
aa,/3(Q;') and the mapping, aa,i3 '■ C{a) — > C{P) is injective. 

Now, suppose that wc have a procedure for computing C'{X), for any given 
P'-closed and bounded semi-algebraic set, X, such that the number and the 
degrees of the polynomials appearing the descriptions of the semi-algebraic 
sets, Xa, a G C'{X), is bounded by 

D''\ (4) 

where Ci > is some absolute constant, and D = J^pev deg(P). 

Using the above procedure for computing C'{X), and the definition of As, we 
have the following quantitative bounds on ^As and the scmi-algcbraic sets 
Xa, a G As, which is crucial in proving the complexity bound of our algorithm. 

Proposition 4.3 Let S (Z he a V-closed semi- algebraic set, where V C 
R[Xi, . . . ,Xjfc] is a family of s polynomials of degree at most d. Then 4^As, 
as well as the number of polynomials used to define the semi- algebraic sets 
Xa, a G As and the the degrees of these polynomials, are all hounded by 
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Proof: Given a e A5 with level(Q;) = j, we first prove by induction on 
level (a) that 

#ancestors(a) < i^nm+-i)-3^j^^)l2 _ 

The claim is clearly true if level(Q;) = 0. Otherwise, from the definition of A5, 
there exists /? G A5, with level(/3) = j — 1, such that a = 70 • • -7^, 7i G 
and m < £ — j + 2. 

For each 7,, we have 

ancestors(7j) =■ ancestors(/9) U {ap^^^i) \ 9 e ancestors(/?)}, 

and it follows that, 

ancestors(Q;) = ancestors(/5) U {a3,6»(7jo) ' ' ' ^i3,e{lin) I 
6 e ancestors (/?), {io, . . . , i„} C {1, . . . , m}}. 

Hence, 

#ancestors(Q;) = # ancestors • 2"* 

< #ancestors(/3) ■ 2^"^'+^ 

< 21^-0(^-^+3) 

< 2^2^', 

for some absolute constant C2. 

We now prove again by induction on the level that there exists an absolute 
constant c > 0, such that the number of elements of A5 of level < j, as well 
as the number of polynomials needed to define the associated semi-algebraic 
sets, and the degrees of these polynomials, are all bounded by {sdY''\ 

The claim is clear for level 0. Now assume that the claim holds for level 
< j. As before, given a e A5 with level(Q;) = j, there exists P E As with 
level(/9) — j — 1, such that a = 70 ■ ■ ■ 7m, 7? £ and m < i—j + 2. We have 
that #ancestors(/3) < 2'^^^ by the previous paragraph. Let ancestors(/3) = 
{9i,...,9n}- Then, 

for 1 < i < iV by the induction hypothesis. 
In order to bound, 

U C'{(3, (3n-(3), 

PieC(ei),i<i<N 
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first observe that < 2^^^^^ and hence the union on the right hand side is over 
an index set of cardinality bounded by, 

and each set in the union has cardinality bounded by, 



where Ci is the constant defined before in (4) above. 

Thus, the total number of admissible indices at level j is bounded by the total 
number of admissible indices at level j — 1 times Z]o<i<€-j+3 (Y) ■ 

It follows that if c chosen large enough with respect to the constants ci,C2, 
then for all k large enough, the the total number of admissible indices at level 
j is at most, 

{sdr. 

The bounds on the number and degrees of polynomials appearing in the de- 
scription can be proved similarly using the same induction scheme. □ 



4-2 Double Complex Associated to a Cover 



Given the different covers described above, we now associate to each a G A5 a 
double complex, and for every /3 e A5, such that a G ancestors 

and leveret) = level a restriction homomorphism: 

satisfying the following: 
(1) 

(Tof(A^''*(a))) = for < i < £ - level(a). (5) 

(2) The restriction homomorphism 

induces the restriction homomorphisms between the cohomology groups: 

for < i < ^ — level(Q;) via the isomorphisms in (5). 
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We now describe the construction of the double complex Ai'''{a) and prove 
that it has the properties stated above. The double complex A4''*{a), is con- 
structed inductively using induction on level(Q;). 

Definition 4.4 The base case is when level(Q;) = i. In this case the double 
complex, M.*''{a) is defined by: 

MP^'^ia) = 0, ifq >0 orp> 1. 

This is shown diagramatically below. 

-0 -0 



-0 -0 



aoSC(a) ao,aiGC(a) 

The only non-trivial homomorphism in the above complex, 
is defined as follows. 

S{(f>)ao,ai = (0ai - 0ao)|x„o.,, for G H°(Xc,J. 

For every jS G A^, such that a G ancestors(/3); and level(Q;) = level(/3) = I, 
we define : M^'^{a) M^'%p), as follows. 

Recall that, M'^'^a) = H°(X„J, and = H°(X^J. 

ao G C(a) Po e C{/3) 

For (j) G M^'^{a) and (3o G C{(3) we define. 
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We define r^'J : M^^^^^a) ^ M^^%p), in a similar manner. More precisely, 
for (f) e M^'°{a) and Po, Pi G C{P), we define 



/3o-/3l ■ 



(The inductive step) In general the A4^''^{a) are defined as follows using in- 
duction on level(Q;) and with ria — i — level(Q;) + 1. 

A^0'9(X) =0, 0<q, 

= ©,„<...<„^, „,ec(a) Tot«(M*'*(ao ■ ■ ■ a^)), < p, < p + q < ria, 
MP''^{a) = 0, else. 

The double complex Ai''* (a) is shown in the following diagram: 



^ Tot"°-^(M*'*(ao -ai)) 
»0«>1 



^ Tot"°-2(X*'*(ao ai aa))- 
aQ<ai<0!2 



^ Tot^(A1*'*(ao-«i)) ^-K ^ Tot=(M'''(ao-ai-a2)) 



^ Tot\M*'*(ao ■ 



a;0<ai<«2 



y 



Totl(M*'*(ao ■ ai ■ c«2) 



3o<ai 



■XQ<ai<<:x2 



H'-fX.^jA 



Tot°(M*'*(ao-c«i)) ^ Tot°(M*''(ao-ai-a2)) ■■■ ^ Tot°(A1*'* (ao • ■ ■ a„„ )) 



D!Q<a;i<a2 



^Q<---<0'no, 



The vertical homomorphisms, d, in Ai*'*{a) are those induced by the differ- 
entials in the various 

Tof{M*''{aQ ■ ■ ■ ap)), ai e C{a). 



The horizontal ones are defined by generalized restriction as follows. Let 
<Pe Tot%M'''iao---ap)), 
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with 

0<j<q 

and 

We define, 

S: Tot'^{M'''{ao---ap)) ^ Tot\M'''{ao ■ ■ ■ a^+i)) 

ao<---<ap,aiGC(a) ao<-<ap+i 



by 



^{4')ao,-,ap+i ( 1) ^ao-cEi-ap+i,ao---ap+i (^ao,.--,cEi,.--,ap+i-^' 

0<i<p+l 0<j<q 



noting that for each i,0 < i < p + 1, ao ■ ■ ■ a, • • • a^+i is an ancestor of 
at) - ■ ■ Op+i, and 



level(Q;o • • • ctj • • • ctp+i) = level(Q;o " " • ctp+i) — level(Q;) + 1, 

and hence the homomorphisms ^^a'o-^i-ap+i.ao -ap+i '^'^^ already defined by in- 
duction. 

Now let, a,(3 & As with a an ancestor of (5 and\eYe\{a) — level(/3). We define 
the restriction homomorphism, 

rXp:M''\a)^M''\P) 

as follows. 

As before, for 4> G M.^'^{a) and (3q & C{(3) we define. 

For < p,0 < p + q < i — level(Q!) + 1, we define 

r^l : ^ M^'^fi), 

as follows. 

Let e M^'^^a) = Tot''(M*''(ao ■ ■ ■ ap)). We define, 

ao<---<ap, ai€C{a) 
(3o<-<l3p,(3iQC{l3) 0<i<g 
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where ap^aiPo ■ ■ ■ Pp) = ai3,a{Po) ■ ■ ■ cii3,a{Pp)- Note that, each ap^aiPi), < i < p 
are all distinct and belong to C{a). Moreover, 

level (a^,c.(^o ■ ■ ■ Pp)) = level(/5o • • • /?p) = level(a) + 1, 

and hence we can assume that the homomorphisms r*'* (^p^...^^-^ i3o-l3p ^^^^ 
the definition of r'J^^ are already defined by induction. 

It is easy to verify by induction on level(Q;) that, M*'*{a) defined as above, 
is indeed a double complex, that is the homomorphisms d and 5 satisfy the 
equations, 

^2 ^ 52 = 0, d 5 + 5 d = 0. 

4-3 Example 

Before proving the main properties of the complexes M.*'*{a) defined above, 
we illustrate their construction by means of a simple example. We take for the 
set S*, the unit sphere S"^ C R'^. Even though this example looks very simple, 
it is actually illustrative of the main topological ideas behind the construction 
of the complex M*'*{S) starting from a cover of S by two closed hemispheres 
meeting at the equator. Since the intersection of the two hemisphere is a topo- 
logical circle which is not contractiblc, Theorem 2.2 is not applicable. Using 
Theorem 2.3 we can compute H°(S'), H^(S'), but it is not enough to compute 
Y{^{S). The recursive construction of M.''* described in the last section over- 
comes this problem and this is illustrated in the example. 




Fig. 2. Example of S"^ C 

Example 4.5 We first fix some notations (see Figure 2). Let Hi and H2 

denote the closed upper and lower hemispheres respectively. Let H12 = Hir\H2 
denote the equator, and let H12 = C1UC2, where Ci, C2 are closed semi-circular 
arcs. Finally, let C12 = Ci fl C2 = {Pi,P2}, where Pi,P2 are two antipodal 
points. 
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For the purpose of this example, we will take for the covers C the obvious 
ones, namely: 



C'{S) 








= {Hi}, i 


= 1,2, 




= {CuC2}, 




C'{Q) 


= {a}, t 


= 1,2, 


C'{Ci2) 


= {Pl,P2}, 






= {Pi}, i 


= 1,2. 



Note that, in order not to complicate notations further, we are using the same 
names for the elements of €'{■), as well as their associated sets. Strictly speak- 
ing, we should have defined, 

C'(S) — {ai, a2}, Xa-^ — Hi, — H2, 

However, since each set occurs at most once, this does not create confusion in 
this example. 

Note that the elements of the sets occurring on the right are all closed, bounded 
contractible subsets of S. It is now easy to check from Definition 4-2, that the 
elements of As in order of their levels as follows. 

(1) Level 

e A5,level(0) = 0, 

and 

C(0) = {ai,a2}, = Hi,Xa2 = H2. 

(2) Level 1: The elements of level 1 are 

ai, a2, ai ■ a2, 

and 

C{a2) = {f32}, Xp, = H2, 

C{ai ■ 02) = {/33,/54}, Xp^ = Ci,Xp^ = €2. 

(3) Level 2: The elements of level 2 are Pi, P2, Ps, Pa, ■ Pa- We also have, 

C(A) = {7i}, X,^ = H„ 1 = 1,2, 
C(A) = {7i}, = Q_2, z = 3,4, 
C(/33 • Pa) = {75, 76}, = Pi-^, ^ = 5,6. 
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We now display diagramatically the various complexes, A4*'*{a) for a e As 
starting at level 2. 



(1) Level 2: For 1 < i < 4, we have 




Notice that for 1 < i < 4, 



The complex A^*'*(/33 • P^) is shown below. 



H°(Pi)0H°(P2 




Notice that, 



H°(Tof (M'''(/33 • m = H°(X^3.^J ^ Q0( 
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(2) Level 1: Fori — 1,2, the complex Ai*'*{ai) is as follows. 









Yi^{Hi) 

Notice that for i = 1,2 and j = 0, 1, 

H^'(Tof(M*''(a,))) = H^'(i/i). 
The complex M.*'*{ai ■ a2)is shown below. 

-0 -0 -0 

-0 -0 -0 

-0 -0 -0 

HO(Ci) 0HO(C2) — H0(P0 0H°(P2) " ^ 

Notice that for j = 0, 1, 



(3) Level 
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The complex Ji4*'*{0) is shown below: 



-0 -0 -0 -0 



-0 -0 -0 -0 



-0 -0 -0 -0 



^H°(Pi)0H°(P2) -0 ^0 ^0 

R\H,)^R\H2) — H°(Ci) 0H°(C2) ^ ^ ^ 

The matrices for the homomorphisms, 5^'*^ and d}'^ in the obvious bases 
are both equal to 

li i\ 

From the fact that the rank of the above matrix isl, it is not too difficult 
to deduce that, ff(Tof(A^''*(0))) ^ ^^{S), for j = 0, 1, 2, that is 

HO(Tof (M''*(0))) ^ Q, 
Hi(Tof(M''*(0))) ^ 0, 
H2(Tof (A4''*(0))) ^ Q. 

We now prove properties (1) and (2) of the various M.*'*{a). 

Proposition 4.6 For each a e As the double complex M.*'*{a) satisfies the 
following properties: 

(1) ff(Tof(A^ •••(«))) ^ ff (X«) /or < i < £ - level(«). 

(2) For every (3 G A^, such that a is an ancestor of {3, and level(Q;) ~ 
level(/3), the homomorphism, r^^'^ : A4''*{a) — > induces the 
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restriction homomorphisms between the cohomology groups: 



r* ■.W{X^)^W{Xp) 

for < i < i — level(Q;) via the isomorphisms in (1). 

The main idea behind the proof of Proposition 4.6 is as follows. We consider 
a triangulation Hq : ^ S, such that for any a E As, ho restricts to a 
semi-algebraic triangulation, ha '■ Aq, Xa- Note that, this implies that if 
(3 E As and a G ancestors(/9), then the triangulation ha '■ Aq, — >• Xa restricts 
to the triangulation /i^ : A^ — > Xf^, and in particular A^ is a subcomplex of 
A«. 

For each a G A^, we have that A^ = Uao£C{a)^aoi ^-iid each A^,, for ctg ^ C{a) 
is a subcomplex of A„. We denote by A/'*'*(A„) the Mayer- Victoris double 
complex of Aq, with respect to the sub-complcxcs Aq,q,q;o G C{a) (cf. Defini- 
tion 3.2). 

Wc denote by = i — level(Q;) -|- 1. Recall that N*''{^a) is the following 
truncated complex. 

Aer'(A„) =Arf''?(A„), 0<p + g<n«, 
M^'^{Aa) = 0, otherwise. 

By Corollary 3.4 we have that, 

(Tof(X:;'(A«))) - W{Xa), 0<i<e- level(a). 

We then prove by induction on level(Q;) that for each a E As there exists a 
double complex D*'*{a) and homomorphisms, 

C : C'(A,) Tof(L>'''(a)) 

such that, 

Tof(0;'*) : Tof{M'''ia)) Tot'{D'''{a)), 

as well as ipa (as shown in the following figure) are quasi-isomorphisms. 
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These quasi-isomorphisms will together imply that, 

ff (Tof(A^*''(a))) ^ (Tof(D'''(a))) = (Tof(AC;'(Ae,))) = (X), 
for < i < £-level(Q;). 

Proof of Proposition 4.6: The proof of the proposition is by induction on 
level(Q;). When levei(Q;) = i, we let D*''{a) = J\f''*{Aa), and define the homo- 
morphisms 0*'*,'0a ^ follows. Prom the definition of A4*'*{a) it is clear that 
in order to define 0*'*, it suffices to define, 0°'° and 0°'^. 

We define, 

ao e C{a) ao £ C(a) 

by defining for 9 e 0^,^^ gc(a) H°(XaQ), and any vertex v of the complex Aq,q, 
^a°(^)ao('^) to be the value of the locally constant function on X^q. 

Similarly, we define 

ao<ai,ai£C{a) ao<aiai£C{a) 

noting that 

ao<ai,ai£C{a) 

and 

A/-r(A„) = C'^(A„„.„J, 

ao<ai,ai eC(a) 

by defining for 9 e 0ao<ai,QieC(a) H'^(Xq,q.q,j), and any vertex v of the complex 
Aq(,.q;i, 0°'^(^)ao,ai ("w) to be the value of the locally constant function 9ao,ai on 
the connected component of X^o-ai containing hao-aiiv)- 

The homomorphism ■0* is induced by restriction as in the definition of ip'^i 
in Corollary 3.4. 

It is now easy to verify, that Tot*(0*'*) and are indeed quasi-ismorphisms. 
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In general for a G A^, with level(Q;) < i, we have by induction that for each 
aQ, .... Op, cip+i G C{a), < p < £ — level (a) + 2, there exists a double complex 
D''*{aQ ■ ■ ■ ap) and quasi-isomorphisms 

Tof(C';...,J : Tof(M'''(ao • • • a,)) Tof(D*''(ao • • • a,)) 
rao-a, : (A„) Tof(L'-'-(ao • • • a,)). 

We now define D*'*{a) by, 

DP''{c^)^eao<-<a,,a,eC{a) Tof^ {D''' {ao ■ ■ ■ a^)) , 0<p + q<na, 
— 0, else. 

The homomorphism 0*'* is the one induced by the different Tot*(0*''..„^) de- 
fined already by induction, that is 

is defined by 

Tot^(C':...aj- 

In order to define the homomorphism ijj^, we first define a homomorphism, 

induced by the different ■0*o - ap- 
We define 

by 

ao<---<cep, aiEC{a) 

We now compose the homomorphism, 

Tof(pr) : Tof(AC:(A„)) ^ Tof(D-'-(a)), 
with the quasi-isomorphism 

C,n„:C;jAJ^Tof(AC:(AJ) 
(see Proposition 3.3). 

Using the induction hypothesis it is easy to see that the homomorphism 
induces an isomorphism between the 'Ei terms of the corresponding spectral 
sequences. It follows from Proposition 3.1 that this implies that Tot*(0*'') is 
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a quasi-ismorphism. A similar argument also shows that Tot'(p''*) is also a 
quasi-isomorphism and hence so is ipa since it is a composition of two quasi- 
isomorphisms. This completes the induction. □ 



5 Algorithmic Preliminaries 



In this section, we describe some algorithmic results which we need in the 
main algorithms. 



5.1 Computation with Complexes 



In the description of our algorithm, we compute in a recursive way certain com- 
plicated double complexes, whose constructions have already been described 
in Section 4. The computation of a complex (or a double complex) means com- 
puting bases for each term of the complex (or double complex), as well as the 
matrices representing the differentials in this bases. Given a complex C* (in 
terms of some fixed bases), we can compute its homology groups II*(C*) using 
elementary algorithms from linear algebra for computing kernels and images 
of vector space homomorphisms. Similarly, given a double complex, D*'*, we 
can compute the complex Tot'{D*'') as well as, H*(Tot*(D*'*)), using linear 
algebraic subroutines. Since the naive algorithms (using say Gaussian elimina- 
tion for computing kernels and images of linear maps) run in time polynomial 
in the dimensions of the vector spaces involved, it is clear that all the above 
computations involving complexes can be done in time polynomial in the sum 
of the dimensions of all terms in the input complex. This is sufficient for 
proving the main result of this paper, and we do not make any attempt to 
perform these computations in an optimal manner using more sophisticated 
algorithms. 



5.2 Ceneral Position and Covers by Contractible Sets 



We first recall some results proved in ()BPR04 ) on constructing singly expo- 
nential sized cover of a given closed semi- algebraic set, by closed, contractible 
semi- algebraic set. We recall the input, ou tput and the complexity of the al- 
gorithms, referring the reader to ( BPR04 ) for all details including the proofs 
of correctness. 
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5.2.1 General Position 



Let Q G R[Xi, ...,Xk] such that Z(Q, R'^) = {x eR'' \ Q{x) = 0} is bounded. 
We say that a finite set of polynomials V C D[Xi, . . . ,Xk] is in strong i- 
general position with respect to Q if any i + 1 polynomials belonging to V 
have no zeros in common with Q in R^, and any i polynomials belonging to 
V have at most a finite number of zeros in common with Q in R*'. 



5. 2. 2 Infinitesimals 

In our algorithms we will use infinitesimal perturbations. In order to do so, 
we will extend the ground field R to, R (e), the real closed field of algebraic 
Puiseux series in e with coefficients in R (|bpro3). The sign of a Puiseux series 



in R(e) agrees with the sign of the coefficient of the lowest degree term in e. 
This induces a unique order on R(e) which makes e infinitesimal: e is positive 
and smaller than any positive element of R. When a G R{e) is bounded by an 
element of R, lim£(a) is the constant term of a, obtained by substituting for e 
in a. We will also denote the field R(£:i) ■ ■ ■ {eg) by R{e), where ei, 62, ■ ■ ■ ,es > 
are infinitesimals with respect to the field R. 



5.2.3 Replacement by closed sets without changing cohomology 

The following algorithm allows us to replace a given semi-algebraic set by a 
new one which is closed and defined by polynomials in general position and 
which has the same homotopy type as the t he give n set. This construction is 
essentially due to Gabrielov and Vorobjov JgvqI), where it was shown that 
the sum of t he Betti n umbers is preserved. The homotopy equivalence property 
is shown in (|BPR04| ). 

Algorithm 5.1 (Cohomology Preserving Modification to Closed) 
Input ; (1) an element c G R, such that c> 0, 

(2) a polynomial Q G R[Xi, . . . ,Xfc] such that Z{Q,R^) C i?(0, 1/c), 

(3) a finite set of s polynomials 

P = {Pi,...,PjcR[Xi,...,X,], 

(4) 0, subset S C Sign{Q,V) , defining a semi- algebraic set X by 

X = U,eE7^(a). 

Output ; A description of a V -closed and bounded semi-algebraic subset, 

X'cZ(g,R(£,£i,...,£2s)'), 

with V = [ji<i<s,i<j<2s{Pi i ^j}) such that, 
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(1) H*(X') = H*(X), and 

(2) the family of polynomials V is in k' -strong general position with respect 
to Z((5,R(e, ^1, . . . ,e2s)^), where k' is the real dimension of 

Z(Q,R(£,£i,...,£2s)'). 

Procedure ; 

Step 1 Let £ he an infinitesimal. 

(1) Define T as the intersection o/ Ext(T, R(£)) with the hall of center 
and radius 1/e. 

(2) Define V as QU {e''{Xl + ... + XI + Xl^^) - 4, X^+i}. 

(3) Replace T hy the V- semi-algehraic set S defined as the intersection of 
the cylinder T x Yi{e) with the upper hemisphere defined hy e'^{Xl + 
... + XI + Xl^,) = 4, Xfc+i > 0. 

Step 2 Using the Gahrielov- Vorohjov construction descrihed in \BPR0^) . re- 



place S hy a V -closed set, S' . Note thatV is in general position with respect 
to the sphere of center and radius 2/e. 

Complexity: Let d be the maximum degree amo ng the p olyuomials in V. 



The total complexity is bounded by (see (|BPR0J)). □ 



5.2.4 Algorithm for Computing Covers hy Contractihle Sets 
The following algorithm described in detail in 

(UproJ is used to a cover 
of a given closed and bounded semi-algebraic sets defined by polynomials in 
general position by closed, bounded and contractible semi-algebraic sets. 

Algorithm 5.2 (Cover by Contractible Sets) 

Input : (1) a polynomial Q eD[Xi, ... ,Xk] such that Z{Q,R^) C B {0,1/ c), 
(2) a finite set of s polynomials V C D[Xi,...,Xfc] in strong i-general 
position on Z((5,R'^). 
Output ; (1) a finite family of polynomialsC = {Qi, ... ,Q]\f} gR.[X I, ... ,Xk] 

(2) the finite family C C R[£][Xi, . . . ,Xfc] (where e denotes the infinitesi- 
mals El ^ 62 ^ ■ ■ ■ ^ e2N > 0) defined hy 

C= {Q±ei\Q EC,l<i< 2N}. 

(3) a set of C-closed formulas {0i, . . . , 0m} such that 
(a) each 7l{(f)i,R{e)'') is contractihle, 

(h) their union Ui<i<Af7?.(0j, R(e)^) = Z{Q,R{e)''), and 
(c) each hasicV-closed suhset ofZ{Q,K{e)'') is a union of some suhset 
of the 7^(0i,R(£)'=)'s. 



Complexity: The total complexity is bounded by s'-^^^^ d'^^''' •* (see (BPRO 
□ 
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6 Algorithm for computing the first i Betti numbers of a semi- 
algebraic set 

We are finally in a position to describe the main algorithm of this paper. 
Algorithm 6.1 [First i Betti Numbers of aV Semi- algebraic Set] 

Input ; a polynomial Q E D[Xi, . . . , Xk] such that 7,{Q, R*^) C B{0,l/c), 

a finite set of polynomials V C D[Xi, . . . , Xk], 

a formula defining a V semi-algebraic set S contained in 7i{Q, R*^). 
Output : bo{S),...,be{S). 
Procedure ; 

Step 1 Using Algorithm 5.1 (Cohomology Preserving Modification to Closed), 
replace S by a V' -closed set, S' . Note that V' is in k' -general position with 
respect toZ{Q,R^). 

Step 2 Use Definition 4-2 to compute As' using Algorithm 5.2 (Cover by 
Contractible Sets) for computing the various C'{-) occuring in the definition 
of As'. For each element a G As', we also compute the set of ancestors 
ancestors(Q;) C As', C{a), as well as level(Q;). 
More precisely, we do the following. 

(1) (a) Initialize, 

As' ^ 0, 

(b) 

As' ^ As' U {0}, 
level(O) ^ 0, 
Xq 5", 

ancestors (0) = {0}. 

Also, maintain a directed graph G with the current set As' as its 
set of vertices representing the ancestor- descendent relationships. 

(2) For i = to i do the following: 

(a) For each a e As' at level i, with ancestors(Q;) = {cti, . . . , a^v}; 

PieC{ai),l<i<N 

using Algorithm 5.2 (Cover by Contractible Sets). 

(b) For 0<j<^ — i + 1 and each ao, . . . , G C{a), 

As' ^ As' U {ao • ai • • • Oj}, 
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■^ao---oij ■* -^ao ri • • • n , 

level(ao ■ ai ■ ■ ■ (Xj) ^ i + I. 

(c) For each {aQ, . . . , ttj} C {/3o, • • • , f^j} C C{a), with j < i — i + I, 

ancestors (/3o ■ ■ ■ Pj) ancestors (/3o • • • U {ao • • • ctj}, 

and update G. 

(d) For each a' e I3n • ex), 

ancestors(Q;') <— ancestors (a') U . . . , /3jv}. 

anc? update G. Use any graph transitive closure algorithm to tran- 
sitively close G. Accordingly update all the sets ancestors(Q;), a e 

Ks>. 

Step 3 Using Definition 4-4) compute for each a G As', the complex Ai''* (a) 
starting with elements a G A^/ with level(Q;) = I. Note that for each a G 
Kg', C{a) has already being computed in Step 2. This allows us to compute 
matrices corresponding to all the homomorphisms in Ai*'*{a) for a G A5/ 
with lcvel(Q;) = £. The recursive definition of A4'''{a) , implies that we can 
compute the matrices corresponding to all the homomorphisms in Ai*''{a) 
for a G As' with level(a) < i, once we have computed the same for M.''* (P) , 
for all (5 G As' with level(/i^) > IcvcI(q). The same is also true for the 
matrices corresponding to the restriction homomorphisms r*'^^. 

Step 4 For each i,0 < i < £, compute 

bi{S) = dimQW{Tof{M'''{0))), 

using standard linear algebra algorithms for computing dimensions of ker- 
nels and images of linear transformations. 

Proof of correctness : The correctness of the algorithm is a consequence 

of the correctness of Algorithms 5.1 (Cohomology Preserving Modification to 
Closed), Algorithm 5.2 (Cover by Contractible Sets), and Proposition 4.6. □ 

Complexity analysis: The complexity of Step 1 is bounded by (sd)*^^'^) us- 
ing the complexity analysis of Algorithm 5.1 (Cohomology Preserving Modi- 
fication to Closed). 

In order to bound the complexity of Step 2, note that the number of calls 
to Algorithm 5.2 (Cover by Contractible Sets), for computing various covers, 
C'(-) is bounded by ^^A^/jWhich in turn is bounded by {sd)k'^^^^ by Proposi- 
tion 4.3. Moreover, the cost of each such call is also bounded by {sd)k'^^^\ 
The cost of all other operations, including updating the list of ancestors of 
elements of A5/ is polynomial in ^As'. Thus, the total complexity of this step 
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is bounded by [sd)k^^^\ Finally, the complexity of the computations involv- 
ing linear algebra in Step 3 is polynomial in the cost of computing the various 
complexes A^*''(a), as well their sizes (see Section 5.1). All these are bounded 
by {sdY°^^^ by Proposition 4.3. Thus, the complexity of the whole algorithm 
is bounded by {sdY ■ n 



7 Implementation and Practical Aspects 



The problem of computing all the Betti numbers of semi-algebraic sets in sin- 
gle exponential time (as well as the related problems of existence of single 
exponential sized triangulations or even stratifications) is considered a very 
important question in quantitative real algebraic geometry. The main result of 
this paper should be considered a partial progress on this theoretical problem. 
Since the complexity of Algorithm 5.2 (Cover by Contractible Sets) for com- 
puting contractible covers is very high (even though single exponential), the 
complexity of Algorithm 6.1 is prohibitively expensive for practical implemen- 
tation. The topological ideas underlying our algorithm has been implemented 
in a very limited setting in order to compu te the first two Betti numbers of 
sets defined by quadratic inequalities (see ( BK05[ )). In this implementation, 
the covering is obtained by means different from Algorithm 5.2. However, 
practical implementation for general semi-algebraic sets remains a formidable 
challenge. 
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